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Abstract 

We present the next-to-leading order QCD corrections to the associated pro- 
duction of the vector gauge boson (Z/W^) and the graviton in the large extra 
dimension model at the LHC. We estimate the impact of the QCD corrections 
on the total cross sections as well as the differential distributions of the gauge 
bosons and find that they are significant. We also study the dependence of 
the cross sections on the arbitrary factorization scale and show the reduction 
in the scale uncertainties at NLO level. Further, we discuss the ultraviolet 
sensitivity of the theoretical predictions. 
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1 Introduction 



The standard model (SM) of particle physics has been very successful in explaining 
the fundamental interactions of the elementary particles, and its predictions have 
been verified experimentally to a very good accuracy except for the discovery of the 
Higgs boson, the only elementary scalar particle in the SM. In spite of its merits, 
the SM has many open questions that are not addressed within its domain and 
a plenty of room is left open for some beyond SM physics scenarios to address 
them. Supersymmetry, extra dimensions, techni-color models are a few to name 
such beyond SM scenarios. With the advent of the high energetic hadron colliders 
like Tevatron and the Large Hadron Collider (LHC), it is quite feasible to probe 
these new scenarios in the laboratory experiments. The LHC with its unprecedented 
center mass energy of 14 TeV and with luminosities as high as 10 34 cm -2 s _1 , offers 
the best possibility of discovering not only the elusive Higgs boson but also of the 
possible new physics that is hidden so far at lower energies. On the other hand, 
the Tevatron which is operating at a center of mass energy of 1.96 TeV has been 
probing such new scenarios. 

One of such beyond SM scenarios that have gained a lot of interest and have been 
studied well in the context of collider phenomenology is the large extra dimensional 
(LED) model proposed by Arkani-Hamed, Dimopoulos and Dvali [lj. This model is 
theoretically well motivated and it addresses the hierarchy problem with the concept 
of extra spatial dimensions. The size of the extra dimensions in this model can be 
of macroscopic size but still consistent with the data from the experiments to date. 
A viable mechanism to hide these extra d spatial dimensions from the SM particles 
is to confine the latter to a 3-brane and allow only the gravity to propagate the full 
4 + d dimensional space time. For simplicity, the extra dimensions can be assumed 
to be flat, of the same size and compactified on a (i-dimensional torus of radius 
Rj (27r). After the compactification, the scale M s of the extra dimensional theory is 
related to the Planck scale M p as: 

M 2 = C d M 2 s +d R d (1) 

where C d = 2 (4tt)-2 /T(d/2) and R is the size of the extra dimensions. This 
compactification implies that a massless graviton propagating in 4 + d dimensions 
manifests itself as a tower of massive graviton modes in 4-dimensions, with mass 
rrii = 4:7r 2 n 2 / R 2 where ft = {ni, n 2 , n d } and = {0,1,2,...}. Here, the zero 
mode corresponds to the 4-dimensional massless graviton. As the inverse square law 
of gravity has been tested down to only few fim so far [2], the size of the extra spatial 
dimensions in this model can be taken as large as this limit. The hierarchy between 
the electroweak scale and the Planck scale can then be accounted for by this large 
volume of the extra dimensions, as can be seen from eqn.([T]). For M s ~ O(TeV), 
the above limit on R constrains the number of extra dimensions to d > 2. 

In the effective theory valid below the scale M s , these gravitons couple to the 
SM fields through energy momentum tensor T^ v of the latter with the coupling 



1 



K = a/167t/M p , as given by [3j H] 



CO 



(2) 



n=0 



Since the coupling is through the energy momentum tensor, the gravitons can couple 
to all the SM fields with the same coupling strength n irrespective of their charge, 
color and flavor. The Feynman rules for the above interaction lagrangian are given in 
[21I1J IS] and in the first reference of [B]. To order k 2 , the above action allows processes 
involving SM fields and virtual gravitons in the intermediate state or real gravitons 
in the final state. In the context of collider phenomenology, this gives rise to a very 
rich and interesting signals that can be seen at the LHC. The virtual exchange of 
the gravitons can lead to the deviations from the SM predictions whereas the real 
emission of the gravitons can lead to the missing energy signals. Though the coupling 
of each graviton mode to the SM fields is M p suppressed, the large multiplicity of 
the available graviton modes can give rise to observable effects. Hence, there will be 
a summation over the graviton modes at the amplitude level for the virtual graviton 
exchanges, and at the cross section level for the real graviton emissions. As the size 
of the extra dimensions could be large in this model, the mass splitting i.e. 2n/R is 
very small and hence this summation over the graviton modes can be approximated 
to be an integral in the continuum limit, with the density of the graviton modes 
given by [I] 



For the real graviton production process at the collider experiments, the inclusive 
cross section is given by the following convolution: 



where da ma is the cross section for the production of a single graviton of mass m^. 
This collective contribution of the graviton modes results in their non-negligible in- 
teraction with the SM fields, and offers the best possibility of probing the low scale 
quantum gravity effects at the colliders experiments. Consequently, a very rich and 
interesting collider signals of some important processes have been reported in the 
literature, but most of them are available only at the leading order (LO) in the 
perturbation theory [3], IU [9] . At the hadron colliders like LHC or Tevatron, 
the QCD radiative corrections are very significant for they can enhance the LO 
predictions as well as decrease the arbitrary scale uncertainties in theoretical pre- 
dictions. Further, the presence of a hard jet in the final state, due to these radiative 
corrections, has the potential to modify the shapes of the transverse momentum 
distributions of the particles that are under study at LO. Obtaining such a modi- 
fication to the shapes of the distributions is beyond the scope of the normalization 
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of the corresponding LO distributions by a constant K-factor, and it requires an 
explicit computation of the cross sections or distributions to next-to-leading order 
(NLO) in QCD. Owing to this importance of the radiative corrections, they have 
been computed for some important processes involving virtual or real graviton ef- 
fects. The K-factors in some cases are found to be as high as a factor of two. Pair 
production processes are the best to exemplify the case of virtual graviton effects, 
where the NLO QCD corrections are computed for di-lepton [6], di-photon [TO], di-Z 
and W + W~ [Uj production processes. In the context of missing energy signals in 
the large extra dimensional model, the NLO QCD corrections are presented for the 
processes (i) jet plus graviton production [13] and (ii) photon plus graviton produc- 
tion [13] . In each of these two cases, it is shown that the K-factors can be as high 
as 1.5 at the LHC. 

In the present work, we compute the NLO QCD corrections to the associated 
production of vector gauge boson and the graviton at the LHC and give a quanti- 
tative estimate of the impact of these radiative corrections. The paper is organized 
as follows. In sec. 2, we discuss the importance of graviton plus vector gauge boson 
production process, outline the phase space slicing method for computing the NLO 
QCD corrections and present the analytical results. In sec. 3, we give the numerical 
results for both the neutral gauge boson and the charged gauge boson cross sections. 
Finally, we present the conclusions in sec. 4. 

2 Graviton plus vector boson production 

The gravitons when produced at the collider experiments escape the experimental 
detection due to their small couplings and negligible decays into SM particles. The 
production of vector bosons (V = Z, W^) together with such an invisible gravitons 
(G) can give rise to a very large missing transverse momentum signals at the collider 
experiments. The study of graviton plus gauge boson production, hence, in general 
will be a useful one in probing the new physics at the LHC. This process has been 
studied at leading order (LO) in the context of lepton colliders [151 HB] as we U as 
at the hadron colliders [17], and also has been implemented in Pythia8 |18] . The 
process is an important one and stands complementary to the more conventional 
ones involving the graviton production, like jet plus graviton or photon plus graviton 
productions, that are generally useful in the search of the extra dimensions at collider 
experiments. 

It is important to note that there is a Standard Model (SM) background which 
gives signals similar to those of associated production of Z and G. This SM back- 
ground receives a dominant contribution coming from the ZZ production process, 
where one of the Z-bosons in the final state decays into a pair of neutrinos (Z — > vv) 
leading to Z-boson plus missing energy signals. The other Z-boson can be identi- 
fied via its decays to leptons, mostly electrons and muons, and then constraining the 
lepton invariant mass close to the mass of the Z-boson to consider only the on-shell 
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Z-bosons. A detailed study of the event selection and the minimization of other SM 
contributions to this process ZZ — > lluu, using MC@NLO and Pythia, is taken up 
in the context of ATLAS detector simulation and is presented in [19J. Any deviation 
from this SM prediction will hint some beyond SM scenario and hence a study of 
this process will be useful in searching the new physics. 

In the context of extra dimensions, a study of the Z plus graviton production at 
LO at the LHC is discussed in [T7j, where the Z-boson identification is done with 
the leptonic decay modes and using the cuts on the leptons as specified in [19J. At 
the LHC, a similar study is done where the signals of Z-boson plus missing energy 
in this model are compared against those coming from the SM [Zvv] background 
and are presented in our recent study [20] . It is worth noting here that a signal of Z- 
boson plus missing energy can also come from the production of Z plus unparticle U, 
where the unparticle leads to missing energy signals. A study of such process based 
on ATLAS detector simulation [17] shows that the vector unparticles are difficult 
to be probed using this channel, whereas the tensor unparticles can give signals 
identical to that of the graviton. In view of the above, it is worth studying the Z 
plus missing energy signals, in particular ZG production, which will be useful to 
confirm the extra dimensional signals once they are seen in the main channels like 
jet or photon plus missing energy. 

In what follows, we describe the computation of NLO cross sections for the 
process under study. Since our focus is on the QCD part in this work, we will 
confine our calculation to the production of on-shell Z-boson. A more detailed study 
involving the Z-boson decays into leptons requires a full detector level simulation 
with the appropriate cuts at NLO and is beyond the scope of the present paper. 



2.1 Leading Order Calculation 

At the lowest order in the perturbation theory, the associated production of the 
vector gauge boson and the graviton takes place via the quark anti-quark initiated 
subprocess, given by 

? (Pi) + Qb(p 2 ) -> V(p 3 ) + G(p 4 ), (5) 

where V — Z, W and a, b are flavor indices. The corresponding Feynman diagrams 
are shown in Fig. ([1]). These diagrams are obtained by considering the tree level 
qqV diagram and by attaching the graviton line to all the possible external legs and 
the qqV vertex. The Feynman rules and the summation of polarization tensor of the 
graviton are given in [3J H]- The couplings of the fermions to the Z and W bosons 
are given by 

-^r(C v -C al 5 ), _^ 7 m(1-7 5 ) (6) 

where 



cos9 w sin9 w ' " y/2 s fn 6 W 
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Figure 1: Feynman diagrams that contribute to the associate production of the 
vector boson and the graviton at the leading order. 

and the coefficients C v and C a are 

C v = Ti-2 sin 2 e w Q f , C a = T( (7) 

Here, Qf and T/ denote the electric charge and the third component of the isospin 
of the quarks respectively, and 9w is the weak mixing angle. For the vector gauge 
boson, the propagator in the unitary gauge (£ — > oo) has been used throughout. This 
choice of the unitary gauge in the electro weak sector has the advantage of having 
vanishing goldstone and ghost contributions. The leading order matrix elements 
for the associated production of Z-boson and the graviton are computed using the 
algebraic manipulation program FORM [21] and are given by (in n-dimensions) 

£W 2 = W^ACl + CDj^- \l2m"(n-2)tu + mHu{3(n-2ft* 



spin 



4 3 96 v v a> (DH 2 u< 

- 2[-68 + n(104 + (-31 + n)n)]t 3 u + 2[284 + n(-264 + (63 - 2n)n)]t 2 u 2 

- 2[-68 + n(104 + (n - 31)n)]tu 3 + 3(n - 2) V - 48m^(n -2)(t + u) 

- 4m^[3(n - 9)(n - 2)t 2 + 2(124 + 3(n - 2l)n)tu + 3(n - 9)(n - 2)u 2 } 

+ 4m2(t + u )[3(n - 5)(n - 2)t 2 - 2(-90 + n(n + 35))tu + 3(n - 5)(n - 2)m 2 ]} 

- 3m 8 {12(n - 2)tit(t + u) + m 2 [(n - 2)V + 2(16 + (n - U)n)tu 

+ (n - 2) 2 u 2 }} - 3m 6 {4m^(18 + (n- 13)n)tu - 2m 2 z {t + u)[(n - 2) 2 t 2 

+ 2(26 + n(2n - 21))tu + (n - 2) 2 u 2 } - tu[(n - 2)(10 + n)t 2 

+ 2(-32 + 3n(2 + n))tu + (n - 2)(10 + n)u 2 }} 

+ 2t 2 u 2 {32m 6 z {n - 2) - 8m^(-6 + n + n 2 )(t + u) 

+ (40 + (n - 17)n) (t + u) [(n - 2)t 2 + 2{n- 4)tu + (n - 2)m 2 ] 
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- m\ \{n - 2) (48 + (n - 25)n)t 2 

+ 2(-156 + ra(118 + (n- 27)n))tu +{n- 2) (48 + (n - 25)n)u 2 ]} 

+ m 4 {48m^(n - 2)tw + 24m 4 (18 + (n - 12)n)tu(t + u) 

- 6tu{t + u) [(n - 2)nt 2 + 2(-12 + n(3n - 4))*u + (n — 2)rax ] 

- m 2 z [3(n - 2) 2 t 4 + 12(n - 6)(3n - 5)£ 3 u + 2(604 + n(25n - 344))tV 

+ 12(n-6)(3n-5)tu 3 + 3(n-2)V]} 

where .D = (s — m|) and s, i, w are the usual Mandelstam invariants. The over all 
bar in LHS of eq. (JSJ) represents that the matrix elements have been averaged over 
the spins and the colors of the initial state particles and summed over those of the 
final state ones. 

2.2 Next-to-Leading Order Calculation 

At the NLO in the perturbation theory, the cross sections receive 0(a s ) contribu- 
tions from real emission as well as virtual diagrams. The integration over the phase 
space of the real emission diagrams will give rise to infra-red (IR) (soft and collinear) 
divergences in the limit where the additional parton at NLO is either soft and/or 
collinear to the initial state partons. On the other hand, the integration over the 
loop momenta in the virtual diagrams will also give rise to infrared divergences, 
in addition to the ultraviolet (UV) divergences. In our calculation, we regulate all 
these divergences using dimensional regularization with n = 4 + e being the number 
of space time dimensions. Completely anti-commuting 75 prescription (22] is used to 
handle 75 in n dimensions. Here, it should be noted that as the gravitons couple to 
the energy momentum tensor of the SM fields, which is a conserved quantity, there 
won't be any UV divergences coming from the loop diagrams. 

There are several methods available in the literature to compute NLO QCD cor- 
rections. Standard methods based on fully analytical computation deal with the 
phase space and loop integrals in n-dimensions and give a finite 0(a s ) contribution 
to the cross sections, after the real and the virtual contributions are added together 
and the initial state collinear singularities are absorbed into the bare parton distri- 
bution functions. However, these methods are not useful whenever the particles in 
the final state are subjected to either histogramming or experimental cuts or some 
isolation algorithms. In such cases, semi analytical methods like phase space slicing 
method or dipole subtraction method are extremely useful. In the present work, we 
have resorted to the former with two cut offs to compute the radiative corrections. 
In this method, the IR divergences appearing in the real diagrams can be handled 
in a convenient way by slicing the soft and collinear divergent regions from the full 
three body phase space. The advantage of this method is that the integration over 
the remaining phase space can be carried out in 4-dimensions, rather than in n- 



dimensions, using standard Monte-Carlo techniques. In what follows, we give some 
of the details about the implementation of this phase space slicing method in our 
NLO computation. 



2.2.1 Real Emission Processes 

There are two types of subprocesses that contribute to the associated production of 
the vector gauge boson and the graviton at NLO in QCD. They proceed by qq and 
qg initial states. At parton level, the 2—7-3 quark anti-quark initiated subprocess 
is given by 

?o(pi) + Qb(p2) -»■ V(p 3 ) + G(p 4 ) + g(p 6 ). 

We find that 14 diagrams contribute to this subprocess and a few of them are 
depicted in Fig. [2J These diagrams are obtained by taking the t-channel qq — > Vg 
diagram at tree level and by attaching the graviton line to all the possible external 
as well as internal lines and to the vertices. The remaining diagrams are obtained by 
interchanging the vector boson and the graviton lines in Fig. [2j In general, diagrams 
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Figure 2: Real gluon emission diagrams 

such as these involving gluons and massless quarks are prone to be singular in the 
soft and collinear regions of the 3-body phase space integration. In the phase space 
slicing method that we adopted here, these soft and collinear regions are separated 
from the full 3-body phase space using two small cut-off parameters, namely S s and 
8 C , that define these singular regions. In the center of mass frame of the partons, 
the soft region is defined as: < E 5 < |<5 s \/s, where E 5 is the gluon energy and y/s 
is the parton center of mass energy. Integration of the eikonal approximated 2 — ?- 3 
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matrix elements over the soft region of the phase space gives the 0(a s ) 2-body 
contribution, 



16 16 

da s = a s C F F(e, fi R , s) ( — + — In 5 S + 8 In 2 5 S ) da (9) 



where 



r(i + |) /4tt^ 



r(i + e) V s 



> = — ^tt — and a s - 



2N 



47T 



Here, a s (/i#) = g 2 (fiR) /47r with p s being the running strong coupling constant, is 
the renormalization scale and N is the number of colors. The region complementary 
to that of the soft "S"', i.e. E 5 > ^S s y/s, is defined as the hard region (H) of the 
phase space. Within this hard region "if", the emitted gluon can be collinear to 
the incoming massless quark or anti-quark and hence can give rise to hard collinear 
divergences. By introducing another small cut-off parameter (<5 C ), we separate these 
collinear divergences from the hard region. The hard collinear region(ifC) can be 
defined as :0 < —tij < S c s (i — 1,2 and j = 5), where £y = (p^ — pj) 2 . In the 
collinear limit, both the 2 — > 3 matrix elements and the 3-body phase space get 
simplified to be expressed in terms of the born cross section as: 



Aa s da F(e,/i R ,s) f ~ jj[P OT (z, e)f q/P (x 1 /z) f q /p(x 2 ) + (q q)] 



dz ( _ 1 - z \ 2 , 
+(x% 4-> X2) } — I o c ) dx\ d%2 



(10) 



where / a /p(^)s' are the bare parton distribution functions (PDF) and P a b(z,e) are 
the unregulated splitting functions in n-dimensions and are related to the usual 
Altarelli-Parisi splitting kernels as P a b(z,e) = P a b(z) — ^,P' ab {z) [23]. Here z denotes 
the fraction of the incoming parton b's momentum carried by the parton a. Note 
that for P qq splitting in the hard region, since a fraction of the parton momentum 
i.e. 8 S is already carried away by the gluon, the effective limits of the integration for 
z will be < z < 1 - 5 S . 

Apart from the qq initiated subprocess at NLO, there will also be a q(q)g initiated 
subprocess given by 



q a {p\) + g{p2) 



V(p 3 ) +G(p 4 ) + q b (p 5 



Here the emitted parton, being a quark or an anti-quark instead of a gluon, won't 
give rise to soft singularity. However, there will be hard collinear singularities when- 
ever the emitted quark (anti-quark) becomes collinear to the incoming partons. 
These collinear singularities are separated using the cut-off 5 C in the same way as in 
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the case of qq initiated subprocess. The cross section in this collinear region turns 
out to be 



dcr HC 9 = Aa s d °0 F ( e , l*R, S ) ( ^j{[ P qg( Z , e ) fq/p( X l) fg/pfa/z) + (? 0)] 

+{x\++xi)j — (s c J dxidx 2 (11) 



These initial state collinear divergences appearing in eqns. (ITOlfe fTTT) as poles in e are 
purely due to the massless nature of the partons involved in the scattering process. 
These divergences can be factored out from the parton level cross sections and be 
absorbed into the bare parton distribution functions at an arbitrary factorization 
scale fip, a process called mass factorization. In the MS scheme, the scale dependent 
parton distribution functions, f a /p(x, fi F ), can be expressed in terms of the bare 
parton distribution functions as given by 



1 

fa/p(x,fi F ) = f a /p{x) + 2a s ^2 Qj F(e,fMR,fMF) J y Pab{z) fb/p( 



x/z), (12) 



where a,b = q(q),g. Substituting these parton densities in da produces collinear 
singular counter terms which when added with the hard collinear contributions 
results in the following 0(a s ) contribution JTUJ ITT] : 

do co u = 2a s da F(e,fJL R ,s)(\f q / p (x2,fi F )[fq/p{xi,fJ>F) + fq/ P {xi, Vf) 

f~ + In -^-J A q^q+g] + (q ?)} + (xi -H- x 2 ] S jdx 1 dx 2 (13) 

where A q ^ q+g = C F (21n5 s + |). The tilde parton distribution functions are given 
by[U[23] 

fq/p{x,fJL F ) = Z2 / —fb/p(x/y,H F )><Pqb(y) (14) 

b=q,g Jx y 



With P ab (y) = P ab (y)\M c l -^-A r )-P' ab {y). (15) 

Note that there is an additional factor of two, as the parton in the final state can be 
collinear to either of the incoming partons, which is implicit from (q <B- q) in eqn. 
f|T3|) . At this stage, one can observe that the divergent pieces that are proportional 
to In 5s cancel among themselves. However, there are singularities still remaining 
that will get cancelled only with those coming from the loop integrals in the virtual 
diagrams. In what follows, we present the details of the virtual corrections to our 
process. 



9 



2.2.2 Virtual Corrections 



The NLO cross sections also receive the contributions coming from the virtual cor- 
rections as well as the wave function renormalization to the 2—^2 leading order 
processes. The corresponding Feynman diagrams are obtained by considering pos- 
sible one loop virtual gluonic corrections to the tree level Feynman diagram for 
qq — > Z and then by attaching the graviton line to all possible internal as well as 
external lines and to vertices, as allowed by the Feynman rules. This way we find 
27 diagrams, out of which 8 diagrams correspond to external leg corrections and 
can be omitted as they vanish in the massless quark limit. Out of the remaining 19 
diagrams, 11 are shown in Fig. [3j The rest of the diagrams can easily be obtained by 
inverting the charge flow direction of the quark lines in the last eight diagrams shown 
in the Fig. ([3]). Interference of these one loop diagrams with the born diagrams gives 




(VWWVWV 



fWWWWV 



Figure 3: Virtual gluon emission diagrams 



0(a s ) contributions. Due to tensorial interaction of gravitons with the SM fields 
the loop integrals involve higher powers of loop momenta in their numerators and 
hence the reduction of tensorial integrals to scalar ones becomes complicated. We 
have written a symbolic program using FORM [21] to perform this reduction in n 
dimensions. The resulting scalar integrals are then evaluated exactly (see [24]) and 
they are listed in the Appendix. Substituting these scalar integrals, we can express 
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the 0(a s ) contribution resulting from the virtual processes as 



16 12 

day = a s da F(e, fi R , s)C F [ + — 



+C 



V x In 2 ( 4 ) + V 2 D f m ( Pl , k, q) + V 3 D> m (P2, k, q) 



fin. 



m 



+V, \n 2 [^)+V 5 In 2 , 



+ V 6 In 2 



-m, 



+ V 7 In 2 



+V 8 In ( -pj + V 9 In + V w i C (k, q) + V n + V 12 + V 13 In j 



+V U In 



m 



+ V 15 (2 



(16) 



where C = cl s k 2 (C 2 + C 2 ) T 2 C F /(AN), C£ n and the D™s are the finite parts of 
the scalar integrals Cq and D respectively and are listed in the Appendix along 
with V^s. It is clear from the above expression that the integration over the loop 
momenta in 4 + e dimensions leads to soft and collinear singularities which appear as 
poles in e. We found that the UV divergences that appear in the intermediate stages 
cancel among various diagrams thanks to the conservation of SM energy momentum 
tensor to this order in perturbation theory. Now, when we add 0(a s ) contributions 
coming from eq. (Q , f|T3l) and ffTBI) , we observe that the remaining soft and collinear 
singularities cancel among themselves as expected, leaving a finite expression for 
the 2-body contribution which can be computed using Monte-Carlo techniques. In 
other words, the 2-body contribution, given by 



da 2 - hody = da, 



da, 



coll 



da 



v- 



(17) 



is found to be free of both UV and IR singularities and hence suitable for further 
numerical evaluation. 

In addition to the above contribution, we also have the hard non-collinear region 
HC of the phase space which do not suffer from any IR singularities by construction. 
The contributions from this region can be obtained by integrating the 2 — > 3 matrix 
elements using standard Monte-Carlo integrations. Owing to the divergence free 
nature of the integration, the 2 — > 3 matrix elements computed in 4-dimensions will 
suffice our purpose. These matrix elements are again computed using FORM. We 
have made several checks to ensure the correctness of our results, namely the gauge 
invariances in QCD, electroweak and gravity sectors. Since contributions from hard 
non-collinear regions involve three body phase space integrals of final state particles 
having different masses, care is needed to parameterize as well as to determine the 
limits of various integrations. We devote our next sub section to discuss this. 
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2.3 Three body contribution 

In this section, we will present briefly how we have implemented various constraints 
imposed by the two cut off phase space slicing method and cuts on the phase space 
integrals for the 2 — > 3 subprocesses. We are interested in the following cross section: 

da 3 - body = J dT,\M 2 q ^ g \ 2 . (18) 

HC,cuts 

where the three-body phase space measure is given by 

dT s = (n^^-)( 27r )^ (4) (Pi+P 2 -P3-P4-P5). (19) 

It is easy to parameterize all the momenta in the center of mass frame of initial state 
partons and then boost them to the lab frame or the center of mass frame of the 
hadrons. The 4-momenta of the massless partons in the initial state, moving along 
the z-axis, are given by 

pi = ^(1,0,0,1), p 2 = ^1(1,0,0,-1) (20) 

where \fs is the parton center of mass energy. The corresponding 4-momenta of 
the massive particles in the final state are given by Pi = (Ei, pi) with masses 
m 2 = E 2 — \pi\ 2 , for i = 3,4,5. For the three body case, it is easy to consider 
the momentum direction of one of the final state particles, say p$, as the reference 
direction and then parameterize the other two momenta p% and p^ with respect to 
this direction: 

p 5 = \p 5 \ {sm6, 0, costf) (21) 

where 9 is the angle between p 5 and the z-axis. The momentum of p% can now be 
parameterized with respect to the direction of p$ and then followed by a rotation in 
the xz-pl&ne by an angle of 9 to get p% = (^3,^3,^3) in the center of mass frame of 
the partons as given by 

P3 — \Ps\ (cos6* cosa sin/3 + sin^ cos/3) 

p\ — [jP3 I sina sin/3 

P3 = \P3\{ cos & cos/3 — sin^ cosa sin/3) (22) 

where a and (3 are the azimuthal and polar angles of p^ with respect to p§. The 
4-momentum of p^ simply follows from the energy momentum conservation. The 
three body phase space in eqn. (fi~9l) can now be expressed in terms of the angular 
variables, using 

^ = d 4 Pl 6(p 2 - m 2 ) = H dEk d^ (23) 
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to get 

dT 3 = dE 3 dtl 3 dE 5 dn 5 5(pl - ml), (24) 

where dfl 3 = <icos/3 da and dQ 5 = dcosd d<j). Further, the angle /3 can be eliminated 
using 

2 \pk\ \M cos/3 = |p 4 | 2 - Wz\ 2 ~ |p 5 | 2 (25) 

Finally, out of the nine integration variables of the three body phase space, in 
eqn. ffl9|) , we are left with four independent variables viz. E 3 , E 5 , 6 and rj, (due 
to 4-momentum conserving delta function and the rotational invariance of P5, the 
reference momentum direction). The three body phase space can then be written 
in terms of these four independent variables as 

dT 3 = g^.^ dE 3 dE 5 dcos9 drj (26) 

The limits of integration of E 3 and E§ can be obtained from the constraint | cos/3 1 < 1, 
and are given by 

E™ = m B , E™ ax = ^=[s + ml-(m 3 + m 4 ) 2 ] (27) 



and 



where 



j-imin,max 

3 " 2r 



A(B + m + m_) ±J(B- m 2 + )(B - m 



(2? 



A = — E5, B = A 2 — \p 5 \ 2 and m± = m 3 ± m 4 (29) 

Finally, all the parton momenta can be boosted back to the lab frame or the center 
of mass frame of the hadrons by a boost factor given, in the limit of the zero rest 
mass of the hadrons, by 

^ = Pern = {Xl-X 2 ) ^ 
E cm (Xl + X 2 ) 

where x\ and x 2 are the fractions of the incoming hadron momenta carried by the 
partons in the center of mass frame of the hadrons. 

We have implemented this phase space parameterization in our numerical code 
written in Fortran 77. We set Pi = ki{i = 3, 4, 5) and 7713 = my, 7714 = m in our code. 
Here my is mass of the gauge boson and m is mass of the graviton. The phase space 
integrations as well as various convolutions in the two and three body contributions 
are done using VEGAS multi dimensional integration package. In what follows, we 
present the impact of our NLO corrections on various observables. 
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3 Numerical Results 



3.1 Neutral gauge bosons 

In this section, we present various kinematic distributions for the associate produc- 
tion of the graviton and the vector gauge boson to NLO in QCD at the LHC The 
results are presented for proton-proton collision energy of y/S = 14 TeV. As dis- 
cussed before, the inclusive cross section for the graviton production involves the 
summation of all possible graviton modes. This summation in the continuum limit 
leads to an integral over the graviton mass. The limits of this integral are set by 
the kinematics from to \/s — my, where \/s is the parton center of mass energy 
and my = mz, my/. The masses of the gauge bosons and the weak mixing angle are 
given by [25] 

m z = 91.1876 GeV, m w = 80.398 GeV, sin 2 w = 0.2312 (31) 

The fine structure constant is taken to be a = 1/128. Throughout our study, we 
have used CTEQ6L1 and CTEQ6.6M parton density sets |26j for LO and NLO cross 
sections respectively. The strong coupling constant is calculated at two loop order 
in the MS scheme with ot s {m z ) = 0.118 (Aq CD = 0.226 GeV). We have also set the 
number of light flavors rif = 5. The following cuts are used for our numerical study, 

p^ w > p™ n , p™ ss > p™ n , \y z ' w \ < 2.5. (32) 

For the 2-body process, the missing transverse momentum is same as that of the 
gauge boson. On the other hand, for the 3-body process, it need not be so due to 
the presence of an observable jet in the final state and hence it amounts purely to 
the graviton transverse momentum. The observable jet is defined as the one that 
satisfies the following conditions: 

p 3 ^ > 20GeV and |r? jet | < 2.5 (33) 

Whenever the jet does not satisfy the above conditions, the missing transverse mo- 
mentum is approximated to be that of the gauge boson. 

The LED model is an effective field theory valid below the UV cut-off scale M s , 
which is expected to be of the order of a few TeV . At the LHC energies 
TeV), it is very well possible that the partonic center of mass energies can exceed 
this scale M s and lead to the signals that do not correspond to the compactified 
extra dimensions of the LED model. This necessitates the need to quantify the UV 
sensitivity of the theory and this issue was already addressed in [3], according to 
which the cross sections can be computed in two different ways, one with truncation 
where the cross sections are set to zero whenever the hard scale Q involved in 
the problem exceeds M s , and the other with un-truncation where there is no such 
constraint imposed on the cross sections. As pointed out in [3j, if these two results 
converge then the predictions are valid and the model is viable, otherwise the un- 
truncated cross sections can dominate the truncated ones, implying the calculations 
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are not under control. In our calculation, we choose the hard scale to be the invariant 
mass of the gauge boson and the graviton, Q = Mzg-, which at LO is the same as 
the center of mass energy of the partons yfs. We will consider both truncated as 
well as un-truncated cases, however, most of our distributions are obtained with our 
default choice of truncation scheme. 

Before proceeding to the kinematic distributions, we will do some consistency 
checks on the calculation. First, we check for the stability of the cross sections 
against the variation of the slicing parameters, 5 S ,5 C . The sum of the 2-body and 
the 3-body contributions given in eqns. ( JT7|) and ffl8|) is expected to be independent 
of the choice of these slicing parameters that are introduced in the intermediate 
stages of the calculation. In fig.fll]), we show the dependency of the transverse 
momentum distribution on the slicing parameter 5 S keeping the ratio of 5 S to 
5 C , fixed at a value of 100. This distribution is obtained using the hard truncation 
scheme for a particular choice of the model parameters M s = 3 TeV and 5 = 4. 
It can be seen from the fig.Q that both the 2-body and the 3-body contributions 
vary with 5 S but their sum is fairly stable against the variation of 5 S over a wide 
range. This ensures the proper implementation of the slicing method in our NLO 
computation. 

Another useful check on the computation is to reproduce the cross sections for 
the associated production of the photon and the graviton at the LHC [14] . In order 
to do this, we recalculated both real emission as well as virtual contributions for this 
process and the corresponding soft and collinear pieces. We found that the following 
replacements 

1. <^±<SUq? 

2. m z 

3. T z -> e 

in the two body and three body real emission matrix elements of the Z/W ± boson 
with Graviton production processes correctly reproduce the corresponding matrix 
elements for photon with Graviton production processes. Here, Qf is the charge 
of the fermion and e is the electromagnetic coupling. Using our symbolic program, 
the analytical expression for virtual contributions for this process agrees with one 
given in the appendix of |l4j . In addition, using these recalculated quantities, we 
reproduced all the numerical results in [H] after taking their choice of parameters, 
cuts etc. It is important to note that the NLO cross sections, or the K-factors, are 
subject to the choice of the event selection or more precisely the cuts on the particles 
in the final state. In our calculation, however, the gauge bosons being massive, we 
present our results according to the cuts given in eqns. (I3"2l fc I3"3l. 

In fig.([n]), the total cross section for the associated production of Z-boson and 
the graviton is shown as a function of P™ n , to NLO in QCD at the LHC. The cross 
sections are given for both the truncated as well as the un-truncated cases and with 
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the choice of model parameters M s = 3 TeV and d — 2. A similar plot is shown 
for d — 4 in fig. OH]). The K-factors are found to have a mild dependency on P™ m , 
varying from 1.6 to 1.4. In fig.([7|), we have shown the variation of the truncated as 
well as un-truncated total cross sections with respect to the scale M s , for the case 
d = 2. The difference between the truncated and the un-truncated cross sections is 
mainly due to the contributions coming from the region Q > M s . However, with 
increasing M s the parton fluxes corresponding to Q in this region rapidly fall down 
and hence the difference between the two cross sections decreases with increase in 
M s . Such a behavior is evident from the figures. ([7]) and ([8]), for d = 2 and d = 4 
respectively. The corresponding K-factors are also shown in fig. ([9]). In the rest of 
our calculation we choose P™ m = 400 GeV and M s = 3 TeV. 

Next, we present LO and NLO transverse momentum distributions of the Z- 
boson (P^) in fig- ffTU]) for d = 2,4,6 and the corresponding missing transverse 
momentum distributions (P™ ss ) in the left panel of fig.ljTT!) for d = 2,4. The QCD 
corrections enhance both P% and P™ ss distributions. Note that the shape of the 
P% distribution remains unaffected while this is not the case for P™ ss distribution. 
Such a pattern can be understood from the definition of P™ lss mentioned before. 

Next, we present the rapidity distributions of Z bosons. The rapidity of massive 
gauge bosons is defined by 



where E and p z are the energy and the longitudinal momentum components of 
the gauge boson in the lab frame. In the right panel of fig. ffTT]) . we have plotted 
the rapidity distribution of the Z-boson both at LO and at NLO for two different 
choices of the factorization scale: Hf = Pt 7 2 and 2P/. This distribution is 
obtained by integrating over the transverse momentum of the Z-boson from 700 
GeV to 750 GeV, for d = 4. Note that the NLO corrections increase the cross 
section. As expected, the inclusion of order a s corrections reduces the dependence 
on the arbitrary factorization scale /i^. The percentage of uncertainty in the cross 
sections at the central rapidity region Y = 0, due the variation of the scale from 
/ip = Pt/2 to fi F = 2P^, is 18.9 at LO and it gets reduced to 8.6 at NLO. 

3.2 Charged gauge bosons 

In this section we discuss the impact of NLO QCD corrections on the associated 
production of charged gauge bosons (Vt^^) and the graviton at the LHC. The matrix 
elements for the W ± case are identical to those for the Z-boson case but for the 
masses of the gauge bosons and their couplings to the quarks as seen in eqn. (J6]). 
Further, in the case of charged gauge bosons, the parton fluxes will also be different 
from those of the neutral gauge boson. The parton fluxes for the quark anti-quark 
annihilation process in the case of Z-boson are of the form qq (q = u, d, s, c, b), while 
they are of the form ud (du) for W + (W~). For W boson production cross sections, 




(34) 
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we will consider the mixing of quarks among different quark generations, as allowed 
by the CKM-matrix elements Vij, with (i = u,c,t) and (j = d,s,b). In view of 
this, in the above parton fluxes, u and d correspond to any -up-type and ciowm-type 
quarks respectively. The CKM matrix elements are given by [25] 

\V ud \ = 0.97425 \V US \ = 0.2252 \V ub \ = 3.89 x 10~ 3 

\V cd \ = 0.230 \V CS \ = 1.023 \V cb \ = 40.6 x 10~ 3 (35) 

Since all our calculations are done in the massless limit of the partons, we have not 
included the top quark contribution in our calculation and set all Vy's to zero. 

Similar to the case of Z boson, we will present the total cross sections as well 
as the differential distributions for the associated production of W boson and a 
graviton. In fig. ffT2]) and fig. ([20]) . we have shown the stability of the transverse 
momentum distributions of W~ and W + respectively, with the slicing parameter 8 S . 
These distributions are obtained for the choice of Pjf = 500 GeV, keeping the ratio 
5 S /5 C fixed at 100. It can be seen from the figures that the sum of the 2-body and 
3-body contributions is fairly stable against the variation of the slicing parameters. 
This ensures the proper implementation of the slicing method in our numerical code, 
taking into account the appropriate parton fluxes for W ± . Next, we present the total 
cross sections as a function of P™" as well as M s . In fig. fll3p and fjg. (jl4j) . we show 
truncated as well as untruncated total cross sections for W~ function of 

P™", for d = 2 and d — 4 respectively. It can be seen from these figures that the 
QCD corrections have enhanced the leading order cross sections considerably, but 
there is no significant change in the shape of the cross sections. Similar plots are 
shown for W + in fig. Q2IJ and fig. Q22p. 

In fig. ffTol) and fig. f|T6l) . we show the total cross sections for W~ as a function 
of M s for d = 2 and d — 4, respectively. A set of similar plots for W + are shown 
in fig. (f23j) and fig. (l24"I) . Note that, in each of the above, the cross sections for W + 
are somewhat higher than the corresponding ones for W~ . This difference in the 
total cross sections can be understood from the respective parton fluxes for W~ and 
W + at the LHC The corresponding K-factors are shown in fig. fll7p for W~ and in 
fig. (|2"5]) for W + . For the choice of the parameters we have considered, the K-factors 
are found to vary from 1.7 to 1.4 in the case of W~ while they range from 1.65 to 
1.05 for W + . Note that the K-factors for W~ case are comparable but a little higher 
than those for W + , which again can be accounted for the differences in the parton 
fluxes. The fact that the valence quark contributions are negligible and the parton 
fluxes at LO for W + are higher compared to those for W~ explains the behavior the 
above factors. 

Further, in figs. (|T8|) and (j26|) . we present the transverse momentum distribution 
of W~ and W + respectively as a function of the number of extra dimensions d and 
for M s = 3 TeV. Similarly, we show the missing transverse momentum distribution 
the graviton when produced in association with W~ in the left panel of fig. (119^1 . In 
the right panel, we present the scale uncertainties in the rapidity distribution of W~ 
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by varying the factorization scale from \ip = P^f /2 to \ip — 2P^ . This rapidity 
distribution is obtained by integrating over the transverse momentum Pjf from 
700 GeV to 750 GeV. Similar plots are shown for W + in fig. fl27]). Note that the 
uncertainty resulting from the factorization scale hf get reduced as we include order 
a s corrections. The percentage of uncertainty at the central rapidity Y w = is 
decreased from 19.1 to 9.3 in the case of W~ , whereas it gets reduced from 18.8 to 
8.3 in the case of W + . 

4 Conclusions 

In this paper, we have systematically computed the full NLO QCD corrections to the 
associated production of the vector gauge boson and the graviton in theories with 
large extra dimensions at the LHC This process plays an important role in probing 
the extra dimensions at the collider experiments, thanks to the large parton fluxes 
available at the LHC. We have used a semi-analytical two cut-off phase space slicing 
method to compute these corrections. We have quantified the ultraviolet sensitivity 
of the theoretical predictions by studying the cross sections in the truncated as well 
as the untruncated cases. In both the cases, the radiative corrections are found to 
have enhanced the cross sections significantly but do not appreciably change their 
shapes. The K-factors for the neutral gauge boson are found to vary from 1.6 to 
1.2 depending on the number of extra dimensions d, while they vary from 1.8 to 1.3 
for the case of charged gauge bosons. Although, the choice of the model parameters 
has the potential to change the cross sections calculated in truncated or untruncated 
cases significantly, we notice that the K-factors remain almost the same in these two 
cases. In addition to the total cross sections, we have also studied the differential 
distributions of the vector gauge bosons and found that the radiative corrections 
are significant and they do not affect their shapes except for the missing transverse 
momentum distribution. At the hadron colliders, the leading order predictions of- 
ten suffer from large uncertainties resulting from the choice of factorization scale. 
Reducing these uncertainties is one of the main motivations for doing NLO compu- 
tation. We have shown that this is indeed the case for the rapidity distributions of 
the gauge bosons by varying the factorization scale from fi F = Pt/2 to /ip = 2Pt, 
leading to reduction in the percentage of scale uncertainty to 9% from 19%. Hence, 
the results presented in this paper are more suitable for studies on associated pro- 
duction of vector boson and graviton in the context of extra dimension searches at 
the hadron colliders. 
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Figure 4: Variation of the transverse momentum distribution of Z boson with 8 S , 
keeping the ratio 5 S /5 C = 100 fixed, for M s = 3 TeV and d — 4. 
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Total cross section for Z and G associated production 
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Figure 5: Total cross section for the associated production of Z and G as a function 
of P^ m at the LHC, for M s = 3 TeV and d = 2. 
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Total cross section for Z and G associated production 
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Figure 6: Total cross section for the associated production of Z and G as a function 
of P^ m at the LHC, for M s = 3 TeV and d = 4. 
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Figure 7: Total cross section for the associated production of Z-boson and the 
graviton at the LHC, shown as a function of M s for d = 2. 
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Total cross section for Z and G associated production 
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Figure 8: Total cross section for the associated production of Z-boson and the 
graviton at the LHC, shown as a function of M s for d — 4. 
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K-factor for total cross section 
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Figure 9: K-factors of the total cross section for the associated production of the 
Z-boson and the graviton at the LHC, given as a function of P™ n (top) and the 
scale M s (bottom). 
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Figure 10: Transverse momentum distribution of the Z-boson for M s 
shown for different values of the number of extra dimensions d. 
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Figure 11: Missing transverse momentum distribution of the graviton produced in 
association with Z-boson at the LHC, for M s = 3 TeV (left). The scale uncertainties 
in the rapidity distribution of Z-boson for M s = 3 TeV and d = 4 (right). 
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Figure 12: Variation of the transverse momentum distribution of W boson with 
5 S , keeping the ratio S S /5 C = 100 fixed, for M s = 3 TeV and d — 4. 



26 



Total cross section for W- and G associated production 
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Figure 13: Total cross section for the associated production of W~ boson and the 
graviton at the LHC, as a function of P™ m for M s = 3 TeV and d = 2. 
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Total cross section for W- and G associated production 
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Figure 14: Total cross section for the associated production of W~ boson and the 
graviton at the LHC, shown as a function of P™ m for M s = 3 TeV and d — 4. 
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Figure 15: Total cross section for the associated production of W boson and the 
graviton at the LHC, given as a function of M s for d = 2. 
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Total cross section for W- and G associated production 
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Figure 16: Total cross section for the associated production of W~ boson and the 
graviton at the LHC, shown as a function of M s for d = 4. 
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Figure 17: K-factors of the total cross section for the associated production of W 
boson and the graviton at the LHC, given as a function of P™ in (top) and the scale 
M s (bottom). 
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Figure 18: Transverse momentum distribution of the W -boson for M s 
shown for different values of the number of extra dimensions d. 



3 TeV is 



10 



10 



10 



10 



Missing transverse momentum distribution 
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Rapidity distribution - Scale variaton 
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Figure 19: Missing transverse momentum distribution of the graviton produced in 
association with ^"-boson at the LHC, for M s = 3 TeV (left). The scale un- 
certainties in the rapidity distribution of W~ boson for M s = 3 TeV and d — 4 
(right). 
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Variation with 8, 
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Figure 20: Variation of the transverse momentum distribution of W + boson with 
5 S , keeping the ratio 5 S /5 C = 100 fixed, for M s = 3 TeV and d — 4. 
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Total cross section for W+ and G associated production 
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Figure 21: Total cross section for the associated production of W + boson and the 
graviton at the LHC, shown as a function of P™ m for M s = 3 TeV and d = 2. 
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Total cross section for W+ and G associated production 
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Figure 22: Total cross section for the associated production of W + boson and the 
graviton at the LHC, shown as a function of P™ m for M s = 3 TeV and d — 4. 
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Total cross section for W+ and G associated production 
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Figure 23: Total cross section for the associated production of W + boson and the 
graviton at the LHC, shown as a function of M s for d = 2. 
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Total cross section for W+ and G associated production 
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Figure 24: Total cross section for the associated production of W + boson and the 
graviton at the LHC, shown as a function of M s for d — 4. 
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K-factor for total cross section 
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Figure 25: K-factors of the total cross section for the associated production of W + 
boson and the graviton at the LHC, given as a function of Prp tn (top) and the scale 
M s (bottom). 
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Figure 26: Transverse momentum distribution of the jy + -boson for M s = 3 TeV is 
shown for different values of the number of extra dimensions d. 
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Figure 27: Missing transverse momentum distribution of the graviton produced in 
association with W^-boson at the LHC, for M s = 3 TeV (left). The scale un- 
certainties in the rapidity distribution of VF+-boson for M s = 3 TeV and d — 4 
(right). 
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Appendix A 

1 Finite part of the virtual contribution 

All the Vi's appearing in eqn. f)16p are given below: 

V 1 = - . \ -((-2m 8 i + 2m 6 (m|(-6t + u) +t(7t + 2u)) -m 4 (l8mU 

(t z u(m z z — s)) 

-2m|(24t 2 + 6tu - u 2 ) + 3t(10t 2 + 8tu + m 2 )) + m 2 t(-8m% + 26t 3 + 36t 2 M 
+I5tu 2 + u 3 + 2m|(21t + hu) - 2m|(30t 2 + 23tw + u 2 )) - At 2 (-2m% + 2t 3 
+4t 2 u + 3tu 2 + u 3 + m%{6t + 4m) - m|(6t 2 + 8tu + 3m 2 )))) 

V 2 = , „ , \ -(I6(2m 8 t + 2m 6 (m 2 (6t-u) -t(7t + 2u)) +m\l8mit 

(t 2 u(m z z — s)) 

-2m|(24t 2 + Qtu - u 2 ) + 3t(10t 2 + 8tu + u 2 )) - m 2 t(-8m% + 26t 3 + 36t 2 u 
+15tu 2 + u 3 + 2m|(21t + hu) - 2m|(30t 2 + 23tu + u 2 )) + 4i 2 (-2m| + 2t 3 
+4t 2 u + 3tu 2 + u 3 + m%(6t + 4m) - m|(6t 2 + 8tu + 3m 2 )))) 

^3 — V2 1 t-^u 
V4 = V\\ t-^u 
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-((AmHu - 2m 6 (9tw(t + u) + m 2 z (t 2 - I2tu + m 2 )) 



(2t 2 u 2 (-m 2 +t + u)) 
+m 4 (36m|tw + 3tw(llt 2 + 16*u + 11m 2 ) + 2m|(t 3 - 30t 2 M - 30tw 2 + m 3 )) 
+4tu(-2m%(t + u) + (t + u) 2 {2t 2 + tu + 2m 2 ) + m|(6t 2 + 8tM + 6m 2 ) 
-m|(6t 3 + llt 2 M + 11£m 2 + 6m 3 )) - m 2 tu(-lQm z + 52m z (t + u) 
-2m|(31t 2 + 46tM + 31m 2 ) + 3(9t 3 + 17t 2 M + lltu 2 + 9m 3 )))) 

(2t¥(-mHH.)) ((4A " 2m6 ( 9tM ^ + M ) + m ^ - l2tu + M ^ 

+m 4 (36m|tM + 3tM(llt 2 + 16tu + 11m 2 ) + 2m|(t 3 - 30t 2 M - 30tM 2 + m 3 )) 

+ttu(-2m%(t + u) + (t + u) 2 {2t 2 + tu + 2m 2 ) + m|(6t 2 + 8tu + 6m 2 ) 

-m|(6t 3 + llt 2 M + lltM 2 + 6m 3 )) - m 2 tu(-16m z + 52m z (t + u) 

-2m|(3H 2 + 46tM + 31m 2 ) + 3(9t 3 + 17t 2 M + 17tM 2 + 9m 3 )))) 
V 6 

((m 2 - t) 2 {m 2 z - t)H 2 {m 2 - t - u)u) ^^(-^z + 
-3m 10 m 2 z (m z (6t - u) + t 2 (Qt + u) - m 2 z t(13t + u)) + m 8 (-27m|t 
+12t 4 M + 3m|(38t 2 + 5tu - u 2 ) + m 2 z t 2 (A8t 2 + lltu + 3m 2 ) 
-m|t(138t 2 + 50tM + 3m 2 )) - 4t 5 M(4m| - 2m|(5i + 4m) 
-t(2t 2 + 5tu + 3m 2 ) + m|(8t 2 + 13tu + 4m 2 )) + m 4 t 2 (12m^° 
-2m|(39t + 7m) + m|(160t 2 + Qltu + 5m 2 ) + t 2 M(56t 2 + 69tM + 12m 2 ) 
+m|t(34t 3 - 107t 2 M - 99tM 2 - 12m 3 ) - 2m|(64t 3 + t 2 M - 2tM 2 + 3m 3 )) 
+m 6 t{-12m 1 z ° - 27m 6 z t(8t + 3m) - 6t 3 M(7t + 4m) + m 8 z (89t + 15m) 
+m|t(-58t 3 + 29t 2 M + 16tM 2 - 2m 3 ) + m|(198t 3 + 94t 2 M + 23£m 2 + 3m 3 )) 
+m 2 t 3 (4m|(4t - u) + 4m z t{8t 2 - 22£m - 17m 2 ) + 4m|(-10t 2 + 6£m + m 2 ) 
-t 2 M(34t 2 + 65tM + 27m 2 ) + 2m|t(-4t 3 + 51t 2 M + 69tM 2 + 18m 3 )))) 

-7—7 5 w \ o 9 7 r^(16(104m 12 mi(t + u) 

(tu(-m 2 + t + u)(-Am 2 m 2 z + (t + u) 2 ) 2 y y zy ' 

+4m 10 m 2 z (96m z (t + m) - 2m|(73t 2 + 74£m + 73m 2 ) - 3(t 3 + 13t 2 M 
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+13tu 2 + m 3 )) + 4(t + uf(-2m z (t + u) + (t + u) 2 (2t 2 + tu + 2m 2 ) 
+m z {6t 2 + 8tu + 6m 2 ) - m|(6t 3 + llt 2 u + Wtu 2 + 6m 3 )) 
+m 4 (3(t + ufillt 2 + 16*it + 11m 2 ) - lQm^iUt 2 + Wtu + 17u 2 ) 
-8m|(t + m) 3 (25£ 2 + 132fu + 25m 2 ) - 12m|(f + m) 2 (43£ 2 - 108tM + 43m 2 ) 
+4m|(t + m) 4 (49£ 2 + 40tw + 49m 2 ) + 48m|(16t 3 + 11£ 2 m + lltw 2 + 16m 3 )) 
-m 2 (t + M) 2 (-16m|(4t 2 + lliu + 4m 2 ) + 2m 2 z (t + u) 3 (9t 2 - 2Uu + 9m 2 ) 
+3(t + M) 4 (9t 2 + 8tM + 9m 2 ) - 4m%(t + u) 2 (A5t 2 + 26tu + 45m 2 ) 
+8m|(25t 3 + 64t 2 M + 64£m 2 + 25m 3 )) + 4m 8 (146m|(t + u) + (t + uf 
+26m|(t + M) 2 (t 2 + 5tu + m 2 ) - 24m|(17t 2 + lQtu + 17m 2 ) 
+4m|(69t 3 + 59t 2 M + 59£m 2 + 69m 3 )) + 2m 6 (136m^°(t + u) - 9(t + m) 6 
-12m|(59t 2 + 54tM + 59m 2 ) - m 2 z (t + m) 3 (123£ 2 + 286tM + 123m 2 ) 
-2m%{t + M) 2 (157t 2 - 374tM + 157m 2 ) + m|(986t 3 + 866t 2 M 
+866tM 2 + 986m 3 )))tt 2 ) 

" (f u \-m^ + t + M)(-4m 2 m| + (t + m) 2 ) 2 ) ((" 96ml2m ^ M 
-At 2 u 2 (t + M) 3 (-6m| + Wm 4 z (t + u) - 15m 2 z (t + m) 2 + 5(t + m) 3 ) 
+4m 10 m|(120m|tM(t + u) + tM<-3t 2 + 34£m - 3m 2 ) + 12m|(t 2 
-\2tu + m 2 )) + m 2 tu{2m 2 z {t + M) 4 (t 2 + 18tu + m 2 ) + (t + u) 5 (3t 2 
+64tM + 3m 2 ) + 2m 4 z (t + u) 3 (9t 2 - 272tu + 9m 2 ) - 4m|(f + M) 2 (llt 2 
-196fu + 11m 2 ) + 24m|(t 3 - 13t 2 M - 13£m 2 + m 3 )) + m A (l92m^tu{t + u) 
-9tu(t + M) 4 (t 2 + 8tM + m 2 ) - 8m 8 z tu(71t 2 + 30tu + 71m 2 ) 
-2m%tu{t + M) 2 (167t 2 - 1060tM + 167m 2 ) + 4m|tM(167t 3 - 227t 2 M 
-227tM 2 + 167m 3 ) - m 2 z (t + m) 3 (3£ 4 - 20t 3 M + 674tV - 20tM 3 + 3m 4 )) 
-2m 8 (432m|tM + tu(t + u) 2 (3t 2 - Atu + 3m 2 ) + 2m|tM(-15t 3 + 179t 2 M 
+ 179tM 2 - 15m 3 ) + 8m|(3t 3 - 97t 2 M - 97£m 2 + 3m 3 ) + 4m|(3t 4 + 118t 3 M 
+66tV + 118tM 3 + 3m 4 )) + m 6 (-384m^°tM + 1312m z tu(t + u) 
+4tu(t + M) 3 (3t 2 + 5tu + 3m 2 ) - 4m|tM(391t 2 + 294tM + 391m 2 ) 
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+m|(f + u) 2 (3t 4 - 70£ 3 m + 1230tV - 70£m 3 + 3m 4 ) + 8m|(3t 5 + 111£ 4 m 

-122tV - 122tV + llltw 4 + 3m 5 )))) 

1 

(3(m 2 - t)t 2 (-m 2 z + t)(m 2 - u)u 2 (-m 2 z + u)(m 2 z - s) 2 (Am 2 m 2 z - (t + m) 2 )) 
x ((-12m 16 m|tM(22m| + 20tu - 1\m\(t + «)) + 3m 14 (20tV(t + m) 2 
+44m|(t 2 - Yltu + m 2 ) - 10m z tu(29t 2 + 118*u + 29m 2 ) + m|(-44t 3 
+852t 2 M + 852tM 2 - 44m 3 ) - 2\m\tu{f - 13t 2 M - Yitu 2 + m 3 )) 
+4tV(t + M) 2 (64m^° - 176m|(t + u) + 8m|(23t 2 + 49tw + 23m 2 ) 
-3tw(7t 3 + 5t 2 M + 5tw 2 + 7m 3 ) - 6m|(16t 3 + 51t 2 M + 51£m 2 + 16m 3 ) 
+m|(24t 4 + lllt 3 M + 134tV + llltM 3 + 24m 4 )) - 3m 12 (792m^°tM 
+4m|(33t 3 - 60K 2 m - 601£m 2 + 33m 3 ) + tV(79t 3 + 257t 2 M + 257£m 2 
+79m 3 ) + m|(-121t 4 + 2026t 3 M + 5262tV + 2026tw 3 - 121m 4 ) 
+m|tM(-73t 4 + 142t 3 M + 574iV + 142tM 3 - 73m 4 ) - m|(llt 5 
+326t 4 M + 3015tV + 3015tV + 326£m 4 + 11m 5 )) - 3m 10 (352m^ 2 tM 
-2428m^°tM(t + u) - t 2 u 2 (t + m) 2 (127£ 2 + 386£m + 127m 2 ) 
-22m|(6t 4 - 167t 3 M - 430tV - 167£m 3 + 6m 4 ) + m|(99t 5 - 1796t 4 M 
-9887fV - 9887fV - 1796£m 4 + 99m 5 ) + m|tM(104t 5 + 521£ 4 m 
+707tV + 707fV + 521tM 4 + 104m 5 ) + m|(33t 6 + 87t 5 M + 2369tV 
+6046t 3 M 3 + 2369tV + 87£m 5 + 33m 6 )) + m 8 (2112m^ 2 tM(t + u) 
-16m^°tM(453t 2 + 1241tM + 453m 2 ) + m|(-132t 5 + 6321t 4 M 
+36839t 3 M 2 + 36839tV + 632Um 4 - 132m 5 ) - 3tV(107t 5 + 809t 4 M 
+ 1912tV + 1912tV + 809tM 4 + 107m 5 ) + m|(33t 6 - 1209£ 5 m 
-20115tV - 42434*V - 20115tV - 1209tM 5 + 33m 6 ) + 3m|tM(84t 6 
+943t 5 M + 3098t 4 M 2 + 4694tV + 3098tV + 943tM 5 + 84m 6 ) + 3m|(33t 7 
-77t 6 M - 544* V + 1792* V + 1792* V - 544* V - 77*m 6 + 33m 7 )) 
+m 2 t 2 M 2 (-1024m^ 2 tM - m|(* + m) 3 (1033* 2 + 3430*m + 1033m 2 ) 
-32m^°(14t 3 - 37* 2 m - 37*m 2 + 14m 3 ) + 3tu(t + m) 2 (41* 4 + 170* 3 m 
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+ 194tV + 170tw 3 + 41m 4 ) + 8m|(142t 4 + 405t 3 M + 478tV + 405tM 3 
+142m 4 ) + m|(483t 6 + 4326t 5 M + 13853t 4 M 2 + 20660fV + 13853£ 2 m 4 
+4326tw 5 + 483m 6 ) - 2m|(69t 7 + 660t 6 M + 2578tV + 5301tV 
+5301tV + 2578tV + 660£m 6 + 69m 7 )) + m 6 (-32m^ 2 to(33t 2 
+131fu + 33m 2 ) + 4m^°tM(537t 3 + 4225t 2 M + 4225£m 2 + 537m 3 ) 
+3*V(t + w) 2 (49t 4 + 469t 3 M + 862tV + 469tw 3 + 49m 4 ) - 4m|tw(126t 4 
+4141t 3 w + 9014tV + 414Km 3 + 126m 4 ) + m|(33t 7 - 630t 6 M + 775tV 
+14222tV + 14222tV + 775tV - 630£m 6 + 33m 7 ) - 3m|tM(43t 7 
+686t 6 M + 3478tV + 7861tV + 7861tV + 3478tV + 686tM 6 + 43m 7 ) 
+m|(-33t 8 + 171t 7 M + 5013t 6 M 2 + 15053t 5 M 3 + 18824tV + 15053fV 
+5013t 2 M 6 + 171tM 7 - 33m 8 )) + m 4 tM(2080m^ 2 tM(t + u) 
-6tu(t + M) 3 (5t 4 + 96t 3 M + 172tV + 96tM 3 + 5m 4 ) + 8m^°(24t 4 
-485t 3 M - 1382* V - 485*m 3 + 24m 4 ) - m|(423* 5 + 375* 4 m - 8054*V 
-8054* V + 375*m 4 + 423m 5 ) + m|(237* 6 + 3942* 5 m + 10563* 4 m 2 
+ 11668fV + 10563* 2 m 4 + 3942*m 5 + 237m 6 ) - m|(39* 7 + 2514* 6 m 
+14408*V + 29959*V + 29959*V + 14408*V + 2514*m 6 + 39m 7 ) 
+m|(33* 8 + 777* 7 m + 5388*V + 16631* 5 m 3 + 24742* 4 m 4 + 16631*V 
+5388tV + 777*m 7 + 33m 8 )))) 

1 

((m| - *) 2 * 2 (m| - u) 2 u 2 {-m 2 + t + u)(-Am 2 m 2 z + (t + m) 2 ) 2 ) 
x (m|(-16m 12 m|tM(6m| - 8m|(* + u) - 2tu(t + u) + m|(3t 2 + 8tu 
+3m 2 )) + 8m 10 m|(6mi?(* 2 - \2tu + m 2 ) + m|(-12* 3 + 149* 2 m + 149*m 2 
-12m 3 ) + 3m|tM(t 3 + 45t 2 M + 45*m 2 + m 3 ) - *V(2* 3 + t 2 u + *m 2 + 2m 3 ) 
+m|tM(3t 4 - 2t 3 M - 54* V - 2*m 3 + 3m 4 ) + m|(6* 4 - 80* 3 m - 284* V 
-80*m 3 + 6m 4 )) + 4*V(* + M) 3 (6m^ 2 - 28m^°(* + u) + 2tu(t + m) 2 (* 2 
+3*m + m 2 ) + m|(49* 2 + HOtM + 49m 2 ) - m|(41* 3 + 155* 2 m + 155*m 2 
+41m 3 ) + m|(17* 4 + 97t 3 M + 166* V + 97*m 3 + 17m 4 ) - m|(3* 5 + 26* 4 m 
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+71tV + 71tV + 26£m 4 + 3m 5 )) - m 8 (864mi 4 tM + 24mi 2 (2t 3 - 117t 2 M 
-117tM 2 + 2m 3 ) - 4m|tw(293t 3 + 1581t 2 M + 1581£m 2 + 293m 3 ) 
+m^°(-72t 4 + 3012t 3 M + 7672fV + 3012tM 3 - 72m 4 ) - 2tV(t 5 + 15t 4 M 
-20t 3 M 2 - 20tV + 15tw 4 + u 5 ) - Am z tu(t 5 + 19t 4 M - 54tV 
-54tV + 19£m 4 + m 5 ) + m 2 z tu(3t 6 + 36t 5 M - 23tV - 688t 3 M 3 
-23tV + 36tw 5 + 3m 6 ) + 2m|(12t 6 + 57t 5 M + 812tV + 1470tV 
+812tV + 57tw 5 + 12m 6 )) - m 2 tM(4m^ 2 (t + M) 2 (23t 2 - 352tM + 23m 2 ) 
-24m^ 4 (t 3 - 13t 2 w - 13tM 2 + m 3 ) + 2tu(t + w) 4 (t 4 + 13t 3 w + 38fV 
+ 13tw 3 + m 4 ) + 2m%(t + w) 2 (59t 4 - 542t 3 w - 2218tV - 542tw 3 + 59m 4 ) 
-2m^°(73t 5 - 845t 4 M - 3284tV - 3284tV - 845tM 4 + 73m 5 ) 
-m|(t + M) 3 (3t 6 + 48t 5 M + 237tV + 236tV + 237tV + 48tM 5 + 3m 6 ) 
+2m z (t + u) 2 (Qt 6 + 87t 5 M + 104t 4 M 2 - 238tV + 104t 2 M 4 + 87tM 5 + 6m 6 ) 
+m|(-49t 7 - 73t 6 M + 2017t 5 M 2 + 6457tV + 6457tV + 2017tV 
-73tM 6 - 49m 7 )) + m 6 (-384m^tM + 193&m%tu(t + u) - 4m^tM(859t 2 
+ 1806tM + 859m 2 ) - t 2 u 2 (t + M) 2 (3t 4 + 64t 3 M - 130tV + 64tM 3 + 3m 4 ) 
+8mf(3t 5 + 347t 4 M + 1095tV + 1095fV + 347tM 4 + 3m 5 ) - m|(45t 6 
+1146t 5 M + 3887t 4 M 2 + 3460t 3 M 3 + 3887tV + 1146tM 5 + 45m 6 ) + m 2 z tu(3t 
+ 115t 6 M + 165tV - 1811t 4 M 3 - 1811t 3 M 4 + 165tV + 115tM 6 + 3m 7 ) 
+2m|(9t 7 + 142t 6 M + 348tV - 1971*V - 197KV + 348tV + 142tM 6 
+9m 7 ) + m|(3t 8 - 32t 7 M - 340t 6 M 2 + 2156tV + 6490£ 4 m 4 + 2156fV 
-340tV - 32tM 7 + 3m 8 )) + m 4 (192m%tu(t + u) - 8m^ 4 tM(101t 2 
+ 138tM + 101m 2 ) + 4m^ 2 tM(333t 3 + 379t 2 M + 379tM 2 + 333m 3 ) 
+2m^°tM(-581t 4 + 30t 3 M + 2118tV + 30tM 3 - 581m 4 ) + 3t 2 u 2 (t + m) 3 (£ 4 
+18t 3 M - 2t 2 u 2 + 18tM 3 + m 4 ) - m 2 z tu(t + M) 2 (3t 6 + 96t 5 M + 209t 4 M 2 
-1132t 3 M 3 + 209tV + 96tM 5 + 3m 6 ) - m|(3t 7 - 613t 6 M + 921tV 
+11961fV + 11961t 3 M 4 + 921tV - 613tM 6 + 3m 7 ) + m|(6t 8 - 177t 7 M 
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+12tV + 8745£ 5 m 3 + 18124tV + 8745tV + 12t 2 M 6 - 177tu 7 + 6m 8 ) 

-m|(3t 9 - 13t 8 u - 338tV + 1474£ 6 m 3 + 8410tV + 8410tV 

+1474tV - 338tV - 13tu 8 + 3m 9 )))) 

1 

((m 2 — t) 2 t 2 (m 2 — u) 2 (m 2 — t — m)m 2 (— 4m 2 m| + (t + u) 2 ) 2 ) 
(m 2 (96m 18 m z tu - Am w m 2 z (158m 2 z tu(t + m) + tu(-3t 2 + 34tu - 3m 2 ) 
+12m|(t 2 - I2tu + u 2 )) - 4*V(* + w) 2 (12m|tM + 6tu(t + m) 3 
+m|(t + w) 2 (t 2 - + m 2 ) - m z {t 3 + 7t 2 M + 7£m 2 + m 3 )) 
+2m 14 (432m|tw + iu(f + m) 2 (3£ 2 - Atu + 3m 2 ) - 4m 2 z tu(8t 3 - 113t 2 u 
-I13tu 2 + 8m 3 ) + 4m|(18t 3 - 281t 2 M - 281£m 2 + 18m 3 ) + 2m|(6t 4 
+411t 3 M + 578fV + 411tw 3 + 6m 4 )) + m 2 t 2 u 2 (A8m 8 z tu(3t 2 - 2tu + 3u 2 ) 
+2(t + w) 4 (3t 4 + 29t 3 M + 72tV + 29tw 3 + 3m 4 ) - m%(t + u) 3 (5t 4 
+124t 3 M - 102tV + 124tM 3 + 5m 4 ) + 2m\{t + m) 2 (7£ 4 + 86t 3 M - 522tV 
+86tM 3 + 7m 4 ) - 4m|(3t 5 + 61t 4 M - 164fV - 164tV + 6Hm 4 + 3m 5 )) 
+m 12 (384m^°tM - 2288m 8 z tu(t + m) - Atu(t + M) 3 (6t 2 + tu + 6m 2 ) - 
16m|(9t 4 - 207t 3 M - 416tV - 207tM 3 + 9m 4 ) - 8m|(9t 5 + 284£ 4 m 
+356t 3 M 2 + 356tV + 284tM 4 + 9m 5 ) - m|(3t 6 - 126t 5 M + 2089tV 
+4756tV + 2089tV - 126tM 5 + 3m 6 )) + m w (-57Qm 1 z °tu(t + u) 
+672m|tM(3t 2 + 7tu + 3m 2 ) + tu(t + u) 2 {39t 4 + 170t 3 M + 222tV 
+ 170tM 3 + 39m 4 ) + 16m|(3t 5 - 147t 4 M - 424tV - 424tV - 147£m 4 
+3m 5 ) + 2m|(36t 6 + 915t 5 M + 780tV - 1286t 3 M 3 + 780t 2 M 4 + 915£m 5 
+36m 6 ) + m|(9t 7 - 109t 6 M + 2141tV + 8903tV + 8903tV + 2141tV 
-109tM 6 + 9m 7 )) - m 4 tM(48m|tM(4t 3 + 3£ 2 m + 3tM 2 + 4m 3 ) 
-m 2 z tu(t + M) 2 (85t 4 - 52t 3 M - 1286tV - 52£m 3 + 85m 4 ) + (f + m) 3 (3£ 6 
+58t 5 M + 303t 4 M 2 + 492t 3 M 3 + 303t 2 M 4 + 58tM 5 + 3m 6 ) - 4m|(3t 6 
+108t 5 M - 55tV - 608t 3 M 3 - 55tV + 108tM 5 + 3m 6 ) + m z {15t 7 
+353t 6 M - 501tV - 5475tV - 5475tV - 501fV + 353tM 6 + 15m 7 )) 
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+m 8 ( y 192m 1 £tu( y t 2 + Atu + m 2 ) - 32m|t«(20t 3 + 97t 2 u + 97tu 2 + 20m 3 ) 
-3tu(t + w) 3 (llt 4 + 80£ 3 m + 122t 2 u 2 + 80tu 3 + 11m 4 ) + 32m|tw(27t 4 
+ 104t 3 M + 77t 2 u 2 + 104tu 3 + 27m 4 ) - 2m|(12t 7 + 423t 6 M + 483tV 
-3970tV - 3970tV + 483fV + 423tw 6 + 12m 7 ) - m|(9t 8 - 36t 7 M 
+880t 6 M 2 + 7376t 5 M 3 + 13062* V + 7376* 3 m 5 + 880* 2 m 6 - 36*m 7 + 9m 8 )) 
+m 6 (-192m^V(t + u) + 16m|*M(3* 4 + 58* 3 m + 62*V + 58*m 3 
+3m 4 ) - 4m|*M(41* 5 + 337* 4 m - 308* 3 m 2 - 308* 2 m 3 + 337*m 4 + 41m 5 ) 
+tu(t + m) 2 (15* 6 + 194* 5 m + 697* 4 m 2 + 1016* 3 m 3 + 697*V + 194*m 5 
+ 15m 6 ) + m|*M(195* 6 + 916* 5 m - 4815* 4 m 2 - 12032* 3 m 3 - 4815*V 
+916*m 5 + 195m 6 ) + m|(3* 9 - * 8 m - 20* 7 m 2 + 2332* 6 m 3 + 7838* 5 m 4 
+7838* V + 2332* 3 m 6 - 20* V - *m 8 + 3m 9 )))) 

Vis = 777^7 q~~r~ m ((144m 10 tM - 6m 8 (97tM(t + m) + 12m|(* 2 - 12*m 

+m 2 )) + 3m 6 (432m|*M + *m(353* 2 + 700*m + 353m 2 ) + m|(48* 3 - 746* 2 m 
-746*m 2 + 48m 3 )) + 12m 4 (48m|*M - 208m|*M(* + u) - *m(89* 3 + 249* 2 m 
+249*m 2 + 89m 3 ) + m|(-6* 4 + 199* 3 m + 482*V + 199*m 3 - 6m 4 )) 
+m 2 tM(-696m|(t + u) + 40m|(39* 2 + 107*m + 39m 2 ) - 6m|(229* 3 
+817* 2 m + 817*m 2 + 229m 3 ) + 21 (27* 4 + 98* 3 m + 126*V + 98*m 3 + 27m 4 )) 
-4*M(-2m|(15* 2 + 86*m + 15m 2 ) + m|(90* 3 + 406* 2 m + 406*m 2 + 90m 3 ) 
-m|(90* 4 + 381* 3 m + 512*V + 381*m 3 + 90m 4 ) + 3(10t 5 + 49* 4 m + 69* 3 m 2 
+69* 2 m 3 + 49*m 4 + 10m 5 )))) (1) 

2 Bq Integrals 



B (V) 



(4tt) 5 



l + 2 - lE -f(V) 



(2) 
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where V G {p3,Pa,P5, k, q}, k is the momentum of vector boson, q, the momentum 
of graviton, p 3 = p 1 + k, p 4 = p 2 + k, p 5 = p 1 + p 2 and, 



f(V) = 




for V = p 3 ,P4 
— in for V = P5, k, q 



3 Cq Integrals 



C (V',V") 



(4tt) 2 [(V - V") 2 - V" 2 } 



-(P' - P") 2 



5 ( T * + lD (£? 



7s + m 



~(p' - v ") 

4:71 flj 



(4) 



where V G {pi,p 2 } and, V" G {A;,g} 



Co(Pi,P2) 



C (k,q) 





% 1 


(An) 2 s 






i\ 


1 6 ~ 




i 1 


(An) 2 s/3 


In 





\Ann. 



-in} + 



2Li 2 

\2 o2 



An^i 2 
2 



27T 



(5) 



- 2Li 7 



1 — a + (3 

'a-p + 1 



In 



a-8-1 



1-a-p 

a + (3 + 1 



In 



a + /3 - 1 



(6) 



where a = and /3 = \^ (t + u) 2 - Am 2 z m 2 

4 D Integrals 

A) (pi,k,q) 



(An) 2 st 
+ (7i? + ln 



4 2 r 



+ ln ( J +ln ( — £ ] +in 

m z z J \m z 



- 7s + In 



47T/i T 2 



+K 7£+hi ( ! 3 +in ( 



- Z7T J - ( 7s + In 
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7T 

~ 12 
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+ -( • > ' I 



3 In 2 



m, 2 — £ 



7T 



-J 



m 



+ 2 In 1 - 



t 



m 2 



+ 1 2in In 1 



t 



2Li, 



+ In 1 



m z 



In 1 - 



In 



77), , 



— 27T 



X 



In 
2 In 
In 2 



m z — s 



+ log 1 - 



m z 



m z 



1 In 1 



m z 



- 2m 



+ ln 2 1 



+ 2 In 



m z 



In 



m 2 



- 1 



(7) 



D (p 2 , k, q) can be readily obtained by replacing l f by V in the above D (pi, k, q) 
expression. 



Do (k,p 2 ,q) 



1 



{An) 2 {tu — m 2 m? 
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-t 



m: 



+ In 



— u 



+ 2ivr 



- 7s + In 
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+ 2Li 2 
+ 2Li 2 
+ In 2 
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2 



— ITT 

An 2 



(m 2 — t) (m 2 — t) 

m 2 m z 2 — tu 

tu — m 2 m 2 
(m 2 — t) (m 2 — u) 

m 2 - t) (m 2 - u)\ 2 / {m 2 - t) {m 2 - u) 



— u 

Airtf 

{m 2 — u) {m 2 — u) 
m 2 m z 2 — tu 
tu — m 2 m 2 
[m 2 — t) {m 2 — u) 



tu — m 2 m z 2 



+ In 



tu — m 2 m z 2 



2 „ f ln (("'-OK -«)) + 1„ ((mS-t)(^-«Y 
\ \ tu — m 2 m z 2 J \ tu — m z m z 2 



(8) 
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Transverse momentum distribution 
n — i — i — i — | — i — i — i — i — | — i — — i — | — i — i — i — i — | — i — i — i — r 
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